Better prediction of protein contact number using a support vector regression analysis of amino acid sequence

Background: Protein tertiary structure can be partly characterized via each amino acid's contact number measuring how residues are spatially arranged. The contact number of a residue in a folded protein is a measure of its exposure to the local environment, and is defined as the number of C-beta atoms in other residues within a sphere around the C-beta atom of the residue of interest. Contact number is partly conserved between protein folds and thus is useful for protein fold and structure prediction. In turn, each residue's contact number can be partially predicted from primary amino acid sequence, assisting tertiary fold analysis from sequence data. In this study, we provide a more accurate contact number prediction method from protein primary sequence. Results: We predict contact number from protein sequence using a novel support vector regression algorithm. Using protein local sequences with multiple sequence alignments (PSI-BLAST profiles), we demonstrate a correlation coefficient between predicted and observed contact numbers of 0.70, which outperforms previously achieved accuracies. Including additional information about sequence weight and amino acid composition further improves prediction accuracies significantly with the correlation coefficient reaching 0.73. If residues are classified as being either contacted or non-contacted, the prediction accuracies are all greater than 77%, regardless of the choice of classification thresholds. Conclusion: The successful application of support vector regression to the prediction of protein contact number reported here, together with previous applications of this approach to the prediction of protein accessible surface area and B-factor profile, suggests that a support vector regression approach may be very useful for determining the structure-function relation between primary sequence and higher order consecutive protein structural and functional properties.

STAR: predicting recombination sites from amino acid sequence

Background: Designing novel proteins with site-directed recombination has enormous prospects. By locating effective recombination sites for swapping sequence parts, the probability that hybrid sequences have the desired properties is increased dramatically. The prohibitive requirements for applying current tools led us to investigate machine learning to assist in finding useful recombination sites from amino acid sequence alone. Results: We present STAR, Site Targeted Amino acid Recombination predictor, which produces a score indicating the structural disruption caused by recombination, for each position in an amino acid sequence. Example predictions contrasted with those of alternative tools, illustrate STAR'S utility to assist in determining useful recombination sites. Overall, the correlation coefficient between the output of the experimentally validated protein design algorithm SCHEMA and the prediction of STAR is very high (0.89). Conclusion: STAR allows the user to explore useful recombination sites in amino acid sequences with unknown structure and unknown evolutionary origin. The predictor service is available from http://pprowler.itee.uq.edu.au/star.

D-amino acid residue in a defensin-like peptide from platypus venom: effect on structure and chromatographic properties

The recent discovery that the natriuretic peptide OvCNPb (Ornithorhynchus venom C-type natriuretic peptide B) from platypus (Ornithorynchus anatinus) venom contains a D-amino acid residue suggested that other D-amino-acid-containing peptides might be present in the venom. In the present study, we show that DLP-2 (defensin-like peptide-2), a 42-amino-acid residue polypeptide in the platypus venom, also contains a D-amino acid residue, D-methionine, at position 2, while DLP-4, which has an identical amino acid sequence, has all amino acids in the L-form. These findings were supported further by the detection of isomerase activity in the platypus gland venom extract that converts DLP-4 into DLP-2. In the light of this new information, the tertiary structure of DLP-2 was recalculated using a new structural template with D-Met(2). The structure of DLP-4 was also determined in order to evaluate the effect of a D-amino acid at position 2 on the structure and possibly to explain the large retention time difference observed for the two molecules in reverse-phase HPLC. The solution structures of the DLP-2 and DLP-4 are very similar to each other and to the earlier reported structure of DLP-2, which assumed that all amino acids were in the L-form. Our results suggest that the incorporation of the D-amino acid at position 2 has minimal effect on the overall fold in solution.

Deduced amino acid sequences surrounding the fusion glycoprotein cleavage site and of the carboxyl-terminus of haemagglutinin-neuraminidase protein of the avirulent thermostable vaccine strain I-2 of Newcastle disease virus

A single-tube RT-PCR technique generated a 387 bp or 300 bp cDNA amplicon covering the F-0 cleavage site or the carboxyl (C)-terminus of the HN gene, respectively, of Newcastle disease virus (NDV) strain 1-2. Sequence analysis was used to deduce the amino acid sequences of the cleavage site of F protein and the C-terminus of HN protein, which were then compared with sequences for other NDV strains. The cleavage site of NDV strain 1-2 had a sequence Motif of (112)RKQGRLIG(119), consistent with an avirulent phenotype. Nucleotide sequencing and deduction of amino acids at the C-terminus of HN revealed that strain 1-2 had a 7-amino-acid extension (VEILKDGVREARSSR). This differs from the virulent viruses that caused outbreaks of Newcastle disease in Australia in the 1930s and 1990s, which have HN extensions of 0 and 9 amino acids, respectively. Amino acid sequence analyses of the F and HN genes of strain 1-2 confirmed its avirulent nature and its Australian origin.

Sulfadoxine resistance in Plasmodium vivax is associated with a specific amino acid in dihydropteroate synthase at the putative sulfadoxine-binding site

Sulfadoxine is predominantly used in combination with pyrimethamine, commonly known as Fansidar, for the treatment of Plasmodium falciparum. This combination is usually less effective against Plasmodium vivax, probably due to the innate refractoriness of parasites to the sulfadoxine component. To investigate this mechanism of resistance by P. vivax to sulfadoxine, we cloned and sequenced the P. vivax dhps (pvdhps) gene. The protein sequence was determined, and three-dimensional homology models of dihydropteroate synthase (DHPS) from P. vivax as well as P. falciparum were created. The docking of sulfadoxine to the two DHPS models allowed us to compare contact residues in the putative sulfadoxine-binding site in both species. The predicted sulfadoxine-binding sites between the species differ by one residue, V585 in P. vivax, equivalent to A613 in P. falciparum. V585 in P. vivax is predicted by energy minimization to cause a reduction in binding of sulfadoxine to DHPS in P. vivax compared to P. falciparum. Sequencing dhps genes from a limited set of geographically different P. vivax isolates revealed that V585 was present in all of the samples, suggesting that V585 may be responsible for innate resistance of P. vivax to sulfadoxine. Additionally, amino acid mutations were observed in some P. vivax isolates in positions known to cause resistance in P. falciparum, suggesting that, as in P. falciparum, these mutations are responsible for acquired increases in resistance of P. vivax to sulfadoxine.

Isolation and characterisation of conomap-Vt, a D-amino acid containing excitatory peptide from the venom of a vermivorous cone snail

Cone snail venom is a rich source of bioactives, in particular small disulfide rich peptides that disrupt synaptic transmission. Here, we report the discovery of conomap-Vt (Conp-Vt), an unusual linear tetradecapeptide isolated from Conus vitulinus venom. The sequence displays no homology to known conopeptides, but displays significant homology to peptides of the MATP (myoactive tetradecapeptide) family, which are important endogenous neuromodulators in molluscs, annelids and insects. Conp-Vt showed potent excitatory activity in several snail isolated tissue preparations. Similar to ACh, repeated doses of Conp-Vt were tachyphylactic. Since nicotinic and muscarinic antagonists failed to block its effect and Conp-Vt desensitised tissue remained responsive to ACh, it appears that Conp-Vt contractions were non-cholinergic in origin. Finally, biochemical studies revealed that Conp-Vt is the first member of the MATP family with a D-amino acid. Interestingly, the isomerization of L-Phe to D-Phe enhanced biological activity, suggesting that this post-translational modified conopeptide may have evolved for prey capture. (c) 2006 Published by Elsevier B.V. on behalf of the Federation of European Biochemical Societies.

Lysine requirements of pre-lay broiler breeder pullets: Determination by indicator amino acid oxidation

The indicator amino acid oxidation (IAAO) method allows the determination of amino acid requirements under conditions of low growth rate as found in pre-laying broiler breeder pullets. Cobb 500 breeder pullets (20 wk old; 2290 +/- 280 g, n = 4) were adapted (6 d) to a pelleted, purified control diet containing all nutrients at greater than or equal to 110% of NRC recommendations. After recovery from surgery for implantation of a jugular catheter, each bird was fed, in random order, test diets containing one of nine levels of lysine (0.48, 0.96, 1.92, 2.88, 3.84, 4.80, 7.68, 9.60 and 14.40 g/kg of diet). Indicator oxidation was determined during 4-h primed (74 kBq/kg body), constant infusions (44 kBq (.) h(-1) (.) kg body(-1)) of L-[1-C-14]phenylalanine. Using the breakpoint of a one-slope broken-line model, the lysine requirement was determined to be 4.88 +/- 0.96 g/kg of diet or 366 +/- 72 mg (.) hen(-1) (.) d(-1) with an upper 95% Cl of 6.40 g/kg of diet or 480 mg (.) hen(-1) (.) d(-1). IAAO allows determination of individual bird amino acid requirements for specific ages and types of birds over short periods of time and enables more accurate broiler breeder pullet diet formulation.

Sequence diversity in the coat protein coding region of the genome RNA of Johnsongrass mosaic virus in Australia

Sequence diversity in the coat protein coding region of Australian strains of Johnsongrass mosaic virus (JGMV) was investigated. Field isolates were sampled during a seven year period from Johnsongrass, sorghum and corn across the northern grain growing region. The 23 isolates were found to have greater than 94% nucleotide and amino acid sequence identity. The Australian isolates and two strains from the U.S.A. had about 90% nucleotide sequence identity and were between 19 and 30% different in the N-terminus of the coat protein. Two amino acid residues were found in the core region of the coat protein in isolates obtained from sorghum having the Krish gene for JGMV resistance that differed from those found in isolates from other hosts which did not have this single dominant resistance gene. These amino acid changes may have been responsible for overcoming the resistance conferred by the Krish gene for JGMV resistance in sorghum. The identification of these variable regions was essential for the development of durable pathogen-derived resistance to JGMV in sorghum.

Endogenous amino acid flow in the avian ileum: quantification using three techniques

The aim of the present study was to compare the protein-free diet, guanidinated casein (GuC) and enzyme hydrolysed casein (EHC) methods for the quantification of endogenous amino acid (AA) flow in the avian ileum. Growing broiler chickens (5 weeks old) were used. All three assay diets were based on dextrose, and in the GuC and EHC diets GuC or EHC were the sole source of N. Endogenous AA flows determined with the use of protein-free diet were considerably lower (P < 0.05) than those determined by the GuC and EHC methods. The, total endogenous AA flows determined by the GuC and EHC methods were almost 3-fold greater (P < 0.05) than those determined by the protein-free diet. The endogenous AA values obtained from GuC and EHC methods were similar (P >0.05), except for the flow of arginine, which was lower (P < 0.05) in the EHC method. Glutamic acid, aspartic acid, threonine and glycine were the predominant endogenous AA present in digesta from the distal ileum. The contents of methionine, histidine and cystine were lower compared with other AA. The method of determination had no effect on the AA composition of endogenous protein, except for threonine, glutamic acid, lysine, arginine and cystine. The concentrations of threonine and arginine were lower (P < 0.05) and that of lysine was higher (P < 0.05) with the EHC method compared with the other two methods. The concentration of glutamic acid was greater (P < 0.05) and that of cystine was lower (P < 0.05) in the EHC and GuC methods compared with the protein-free diet method. The results showed that the ileal endogenous flows of N and AA are markedly enhanced by the presence of protein and peptides, above those determined following feeding of a protein-free diet. It is concluded that the use of EHC and GuC methods enables the measurement of ileal endogenous losses in chickens under normal physiological conditions.

Sequence, circulating levels, and expression of C-type natriuretic peptide in a euryhaline elasmobranch, Carcharhinus leucas

The present study has examined expression and circulating levels of C-type natriuretic peptide (CNP) in the euryhaline bull shark, Carcharhinus leucas. Complementary DNA and deduced amino acid sequence for CNP in C leucas were determined by RACE methods. Homology of CNP amino acid sequence in C. leucas was high both for proCNP and for mature CNP when compared with previously identified elasmobranch CNPs. Mature CNP sequence in C. leucas was identical to that in Triakis seyllia and Seyliorhinus canicula. Levels of expression of CNP mRNA were significantly decreased in the atrium but did not change in either the brain or ventricle following acclimation to a SW environment. However, circulating levels of CNP significantly increased from 86.0 +/- 7.9 fmol ml(-1) in FW to 144.9 +/- 19.5 fmol ml(-1) in SW. The results presented demonstrate that changes in environmental salinity influences both synthesis of CNP from the heart and also circulating levels in C. leucas. Potential stimulus for release and modes of action are discussed. (c) 2005 Elsevier Inc. All rights reserved.

The crystal structure of a bacterial Class II ketol-acid reductolsomerase: Domain conservation and evolution

Ketol-acid reductoisomerase (KARI; EC 1.1.1.86) catalyzes two steps in the biosynthesis of branched-chain amino acids. Amino acid sequence comparisons across species reveal that there are two types of this enzyme: a short form (Class 1) found in fungi and most bacteria, and a long form (Class 11) typical of plants. Crystal structures of each have been reported previously. However, some bacteria such as Escherichia coli possess a long form, where the amino acid sequence differs appreciably from that found in plants. Here, we report the crystal structure of the E. coli enzyme at 2.6 A resolution, the first three-dimensional structure of any bacterial Class 11 KARI. The enzyme consists of two domains, one with mixed alpha/beta structure, which is similar to that found in other pyridine nucleotide-dependent dehydrogenases. The second domain is mainly alpha-helical and shows strong evidence of internal duplication. Comparison of the active sites between KARI of E. coli, Pseudomonas aeruginosa, and spinach shows that most residues occupy conserved positions in the active site. E. coli KARI was crystallized as a tetramer, the likely biologically active unit. This contrasts with P. aeruginosa KARI, which forms a dodecamer, and spinach KARI, a dimer. In the E. coli KARI tetramer, a novel subunit-to-subunit interacting surface is formed by a symmetrical pair of bulbous protrusions.

Discovery of an MIT-like atracotoxin family: Spider venom peptides that share sequence homology but not pharmacological properties with AVIT family proteins

This project identified a novel family of six 66-68 residue peptides from the venom of two Australian funnel-web spiders, Hadronyche sp. 20 and H. infensa: Orchid Beach (Hexathelidae: Atracinae), that appear to undergo N- and/or C-terminal post-translational modifications and conform to an ancestral protein fold. These peptides all show significant amino acid sequence homology to atracotoxin-Hvf17 (ACTX-Hvf17), a non-toxic peptide isolated from the venom of H. versuta, and a variety of AVIT family proteins including mamba intestinal toxin 1 (MIT1) and its mammalian and piscine orthologs prokineticin 1 (PK1) and prokineticin 2 PK2). These AVIT family proteins target prokineticin receptors involved in the sensitization of nociceptors and gastrointestinal smooth muscle activation. Given their sequence homology to MITI, we have named these spider venom peptides the MIT-like atracotoxin (ACTX) family. Using isolated rat stomach fundus or guinea-pia ileum organ bath preparations we have shown that the prototypical ACTX-Hvf17, at concentrations up to 1 mu M, did not stimulate smooth muscle contractility, nor did it inhibit contractions induced by human PK1 (hPK1). The peptide also lacked activity on other isolated smooth muscle preparations including rat aorta. Furthermore, a FLIPR Ca2+ flux assay using HEK293 cells expressing prokineticin receptors showed that ACTX-Hvf17 fails to activate or block hPK1 or hPK2 receptors. Therefore, while the MIT-like ACTX family appears to adopt the ancestral disulfide-directed beta-hairpin protein fold of MIT1, a motif believed to be shared by other AVIT family peptides, variations in the amino acid sequence and surface charge result in a loss of activity on prokineticin receptors. (c) 2005 Elsevier Inc. All rights reserved.